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ABSTRACT 

We apply the stochastic model of iron transport developed by Rebusc o" et alJ 1*2005) to the 
Centaurus cluster. Using this model, we find that an effective diffusion coefficient D in the 
range 2 x 10 28 — 4 x 10 28 cm 2 s _1 can approximately reproduce the observed abundance dis- 
tribution. Reproducing the flat central profile and sharp drop around 30 — 70kpc, however, 
requires a diffusion coefficient that drops rapidly with radius so that D> 4 x 10 28 cm 2 s _1 
only inside about 25kpc. Assuming that all transport is due to fully-developed turbulence, 
which is also responsible for offsetting cooling in the cluster core, we calculate the length 
and velocity scales of energy injection. These length scales are found to be up to a factor of 
~ 10 larger than expected if the turbulence is due to the inflation and rising of a bubble. We 
also calculate the turbulent thermal conductivity and find it is unlikely to be significant in 
preventing cooling. 
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1 INTRODUCTION 

Galaxy clusters contain bright, extended, sources of X-ray emis- 
sion. This emission is primarily thermal bremsstrahlung in the low- 
density (~ 10~ 3 cm~ 3 ), optically-thin plasma of the IntraCluster 
Medium (ICM), gravitationally heated to 10 7 — 10 8 K in the clus- 
ter dark matter potential. Spectroscopic observations have revealed 
the ICM to be enriched with metals ejected from the cluster galax- 
ies with an avera ge background metallicity of about a third solar 
I Edge ~Stewartfll99ll) . Observ ations show that this value is con- 
stant out to a redshift of z ~ 1 I Mus hotzkv & Loewensteinlll997t 
lHashimoto et al J [20041) and a comparison of the abundance dis- 
tribution with supernova models has led to the conclusion that 
this background enrichment is the result of type II supernova ac- 
tivity early in th e cluster's history (e.g. iFinoguenov et alJl200Ct 
iTozzi et all2003l) . 

Two distinct classes of cluster have been identified based on 
their X-ray surface brightness profile. Cool-core clusters are char- 
acterised by a bright peak in their X-ray surface brightness profile, 
which is not present in non-cool-core clusters. The cool-core clus- 
ters are believed to correspond to dynamically relaxed systems, i.e. 
those which have not undergone a recent merger. The ICM cooling 
timescale in these bright central regions is often much shorter than 
the Hubble time. It was expected that this short cooling time would 
result in the deposition of large amounts of cool gas that would be 
obser ved in these regions — the cooling flow model (e.g. [Fabian 
1994). However, detailed observations of clusters have failed to 
locate gas below about a third of the mean cluster temperature 
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(e.g. lAllen et alJl200ll: [Peterson et aill2003l Islanders et aljEool 
IVoigt & Fabianl2004 and it is now widely believed that some form 
of distributed heating must be taking place. 

A wide range of solu tions to the cooli ng flow problem have 
been proposed (see lPeterson & F abian 2005, for a review), with 
most of the attention focused on energy injection by the central 
Active Galactic Nucleus (AGN). These models have the advantage 
that the AGN is able to prov ide enoug h power to offset the cooling 
flabor & Binne^ll993t Ichurazov et alJl2002t : however, it is un- 
clear how the energy can be distributed to provide enough heating 
over the entire cooling region. A clue may be provided by the ob- 
servation of bubbl es of sync hrotron-emitting plasma i n many cool- 
core clusters (e.g. |BJjhrjngex^LalJLL993|; [Sanders & Fabianll2002t 
iMcNamara et alJl200Ct iBirzan et alj|2004 . These are inflated by 
jets emanating from the central AGN. Some cl usters show ev idence 
for further "ghost" bubbles at larger radii (e.g. lFabian et all2002l) . 
These show weak, low-frequency, radio emission and are presum- 
ably the result of an earlier epoch of AGN activity. This has led to 
much work on the possibility that the inflation and subsequent rise 
of these bubbles through the cluster atmosphere is responsible for 



the transfer of eners 


y from AGN to ICM (e.s 


iMathews et all200^ 


.iRuszkowski et alj2004) 



There is speculation that cluster heating might be driven by 
turbulent motion of the cluster gas. Evidence for a Kolmogorov- 
type turbulent casc ade has b een found in the pressure fluctuation 
spectrum of Coma I Schuecker et al. 2004) and in the Faraday rota- 
tion power spectrum of Hydra A l Vogt & EnBlin 2005). However, 
the level of turbulence in cool-core clusters remains unclear; the 
observation of extended, linear Ha filaments in several cool-core 
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clusters has been used as evidence that the ICM in these systems 
cannot be fully turbulent ( Fabia n et all2003bl) . 

A method of quantifying the level of stochastic motions in 
cool-c ore clusters wa s devel oped and applied to the Perseus clus- 
ter bv iRebusco et alj (2005). Their method is based on the cen- 
tral metallicity peak observ ed in cool-core clusters which is absent 
in no n-cool-core systems iFukazawa et aljEoo4 iDe Grandi et alJ 
2004). This is likely a result of metal ejection by the Bright- 
est Cluster Galaxy (BCG). The abundance profile does peak on 
the BCG, but does not match the light profile in detail. By as- 
suming the difference in these profiles to be the result of diffu- 
sion of metals cause d by stochastic motions in the cluster gas, 
IRebusco et alJ 120051) were able to determine the diffusion coeffi- 
cient; D = 2 x 10 29 cm 2 s _I for the Perseus cluster. Linking this 
with a model of turbulent heating, and assuming heating balanced 
cooling, they predicted the length and velocity scales of the gas 
motion as / ~ lOkpc and v ~ 300 km s . 

Here we apply the model of lRebusco et alJ J2005) to the case 
of the Centaurus cluster. Centaurus is unusual in showing super- 
solar abundances near the centre - the peak metallicity is around 
2 compared with about 0.7 for Perseus - making it an interesting 
limiting case for the model. 

We adopt the cosmology Hq = 70kms~ 1 Mpc -1 , Q. m = 0.3, 
Qa = 0.7 and use the solar abundances of lAnders & Grevessd 
1 1989) throughout. 



2 DIFFUSION MODEL 

We employed the model for iron transport developed in 
IRebusco et alJ 120051) . This is based on a combined diffusion- 
advection equation in which the advection term ensures that the 
equations satisfy the physical constraint that a spatially uniform 
abundance undergoes no evolution in the absence of sources: 

^^=V-(n H DV(d))+S, (1) 

tin is the ICM hydrogen density, a is the abundance, D is the diffu- 
sion coefficient and 5 the (spatially distributed) so urces of iron. 

The model of iron injection is taken from iBohringer et alJ 

(2004) and considers two sources of iron; type la supernovae and 
stellar winds. In each case, it is assumed that the rate of injection 
was greater in the past when the stellar population of the BCG was 
younger and that the chang e in time can b e approximated using the 
power law model of iRenzini et alJ ^93) leading to source terms 
of the form: 

S SN = 0.7xl0- I2 xStfx (J-J (^d^M Q yr-\ (2) 

SsteUar = 1 X ^ (^) Mq^K (3) 

Here L b q and Mq refer to the solar B-band luminosity and 
mass, respectively, f/, is the Hubble time, SR is the supernova rate 
in SNU (one SNU is the rate of supemovae explosions such that a 
galaxy with 10 10 L b q has one SN per century), k controls the rate 
of type IA supernovae in the past and Lb describes the luminos- 
ity profile of the galaxy. The model parameters k, SR and the total 
integr ation time are no t known. We expect k in the range 1.1—2 
( Bohri nger et aljEo04l) . The integration time is interpreted as the 



time since the last major cluster merger, as such a process would 
have a disruptive effect on the abundance peak. This time together 
with the values of k and SR are subject to the constraint that the 
total mass injected match the observed excess in the cluster. Fol- 
lowing Rebusco et al. (2005), we take k = 1.4, integrate for 8Gyr 
and fix the supernova rate at 0. 19 SNU to match the observed mass 
excess. 

In addition to iron injection, there will be evolution of the ICM 
hydrogen density with time; for example, stellar mass-loss liberates 
hydrogen as well as heavy elements, and merging subsystems will 
have gas stripped. We neglect these effects and regard the hydrogen 
density as a constant in all our models. 



3 MODELLING DIFFUSION IN CENTAURUS 
3.1 Properties of the Centaurus Cluster 

High resolution X-ray imaging of the core of the Centaurus clus- 
ter has revealed a complex, asymmetric, structure with a wealth of 
small-scale features (Fabia rTet all2005l) . For our purposes it is suf- 
ficient to consider a spherically symmetric profile. In order to con- 
struct an accurate profile over a wide range of cluster radii, we have 
combined data from several observations of Centaurus. The high- 
est resolution Chandra data is only available for the inner ~ 60kpc 
of the cluster and so we have supplemented this with a new reduc- 
tion of archival XMM-Newton data out to ~ lOOkpc and published 
deproj ections of d ata from the ROSAT (Allen & Fabian 1994J) and 
ASCA <Whitel200d) satellites at larger radii. 

3.1.1 Data Reduction 

We analysed Chandra observation IDs 0504, 0505, 4954, 4955 
and 5310. Each observation placed the centre of the cluster near 
the aim-point of the ACIS-S3 CCD. The dataset was reprocessed 
using CIAO and the latest appropriate gain file (acisD2000-01- 
29gain_ctiN0003.fits). Flares in the datasets were filtered by using 
the 2.5 - 7keV count rate on the ACIS-S1 CCD (which is back- 
illuminated like the S3 CCD), and the LCLCLEAN tool. This yielded 
a total exposure time of 196 ks. The event files were reprojected to 
match the longest observation (4954). Standard blank-sky observa- 
tions were tailored to match the individual observations, and then 
reprojected to match 4954. The exposure time in each background 
observation was adjusted so that the count rate in the 9 — 12keV 
band matched the cluster observation, to account for the gradual 
increase in background with time. 

In each annulus, spectra and background spectra were ex- 
tracted from each of the datasets and blank-sky observations. 
Weighted response matrices were generated for each of the dif- 
ferent observations using MKACISRMF, and weighted ancillary re- 
sponses created using MKWARF. The extracted spectra in the annu- 
lus were added together. We took the background spectra for each 
of the datasets. The exposure time of the spectra were effectively 
reduced by discarding photons. The exposure times were reduced 
so that the fractional time of each background observation to the 
total exposure time of the background spectrum matched the frac- 
tional time of each observation to the total exposure time of the total 
foreground spectrum. This procedure also optimised for the longest 
total background exposure time. The shortened background spectra 
were added together to create a total background. The responses 
were added together, with weighting according to the number of 
counts as before. 
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The spectra were fit using the PROJCT mod el in XSPEclAr naud 

ll996h to account for projection. A VMEKAL iMewe et alJll985t 
Lieda hl et alll99.^) model was used to compute the emission spec- 
tra. In the fit the Galactic absorption was a free parameter. In each 
shell, the temperature, normalisation, O, Ne, Si, S, Ar, Ca, Fe and 
Ni abundances were free parameters. The model was fit between 
0.5 and 7keV, minimising the % 2 statistic. 

To obtain the ID profile from the Chandra observations, we 
averaged over a 93.9-degree sector to the west of t he centre, s o 
avoiding the swirl of cool gas seen to the east l Fabian et all 2005). 

The XMM-Newton data were reduced in standard fashion. 
High-energy (10- 15keV), single pixel (PATTERN == 0) light 
curves were constructed in 100-second bins for the events satis- 
fying the FLAGs #XMMEA_EM / #XMMEA_EP for the MOS / pn, re- 
spectively. Count-rate cuts of 0.1 5/0.45 ct s _1 were applied for the 
two instruments. This results in 3 1 ks of good time for the MOS, 
but zero for the pn due to background flaring. 

The region between — 150kpc was divided into 10 annuli 
with equal net counts. Backgroun d spectra were extracted from the 
blank-sky fields of Read & Ponman ( 2003), and normalised using 
the ratio of the count-rate outside the field of view for the science 
to background observations. The results presented are from MOS 1 
only. The model was pro j ct*phabs (mekal) with a common, free 
abundance, and spectra were fitted over the range 0.5 — 7.5keV. 



3.1.2 Adopted profiles 

The deprojected iron abundance profile used is shown in Fig.Q We 
have removed data-points corresponding to a central abundance de- 
crease. Such an "abundance hole" could be due to the gas being 
multiphase. If it is real, it will certainly require some additional 
physical process beyond diffusion from a centrally concentrated 
sourc e, e.g. the radial gas motion associated with bubble forma- 
tion IChurazov et alj |2004). We have fitted a model profile of the 
form: 



19.2+ te) 4 ' 55 

a = 0.3 \ / - T rrap ) . 

2.95 + (^) 4 - 55 



(4) 
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Figure 1. Deprojected iron abundance profile for the Centaurus cluster. 
The data is from Chandra (circles), XMM-Newton (diamonds) and ROSAT 
(squares). Chandra data for the inner 5 kpc (where an abundance hole is 
observed) has been excluded. The solid line represents a fitted profile with 
a background abundance of 0.3 solar, adopted through the remainder of the 
work. 



The deprojected hydrogen density profile is shown in Fig. [2] 
with the adopted fit of the form: 
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Figure 2. Deprojected hydrogen density with data from Chandra (circles), 
(5) XMM-Newton (diamonds) and ROSAT (squares). The solid line is a power- 

law fit 



It has been assumed that the electron and hydrogen densities are 
related as n e = 1 .Inn- The deprojected temperature profile is shown 
in Fig.|3| with a powerlaw fit of the form: 



/ \ 0.38 

T = 0.72 [ — ) keV, 
which is valid over the range 1 < r < 100 kpc. 



for absorption, and B = 1 1.03 with the stated correction for galac- 
tic B-band extinction. Hence we adopt a B-band luminosity for 
(6) NG C 4696 of 1 .2 x 1 n L w . The luminosity profile was taken to 
be a Hernauist ( 1990) profile with a scale radius of 9.9kpc chosen 
to match the radius c ontaining 50 per cent of the total B light in 
Lauberts & Valentiidil989l) . 



3.2 Properties of NGC 4696 

NGC4696 is the brightest cluster galaxy of the Centaurus clus- 
ter and is situated at the centre of the abundance peak. We derive 
the luminosity of NGC 4696 from the integrated photometry pre- 
sented in the HYPERLEDA catalogue iPaturel et all2003l) . The to- 
tal B-band magnitude is given as B = 11.51 ±0.01, uncorrected 



4 RESULTS 

Fig. [5] shows the abundance profile produced by simulations us- 
ing spatially uniform diffusion coefficients of 2 x 10 28 cm 2 s _1 , 
4 x 10 28 cm 2 s _1 and 2 x 10 29 cm 2 s~'. For comparison, Fig. |4] 
shows the abundance profile in the absence of diffusion. D = 
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Figure 3. Deprojected temperature profile for the Centaurus cluster with 
data from Chandra (circles), XMM-Newton (diamonds) and ASCA (stars) 
The solid line is a simple power-law fit to the profile. This form provides a 
good fit in the inner 100 kpc of the cluster but may not be accurate outside 
this region. 




Radius / kpc 



Figure 4. Background-subtracted abundance profile in the absence of any 
diffusion; i.e. when the abundance follows the BCG light. Points represent 
the data, as in figureQ] whilst lines represent the model at times of 1 ,2,4 and 
8 Gyr. 



4 x 10 28 cm 2 s _1 provides a reasonable fit to the observed abun- 
dance profile in the central regions, although the model signifi- 
cantly underestimates the iron abundance at higher radii, whilst 
D — 2 x 10 28 cm 2 s provides a better fit at outer radii but is a 
poor fit within 20 kpc of the core. These values of D are factors 5 



Perseus by Reb usco et alj 120051) and are consistent with the limit 
D < 6 x 10 28 cm 2 r' derived bv lFabian etal] l200g) based on the 
width of the abundance peak and the local sound crossing time. 



4.1 Variable diffusion coefficient 

iRebusco et aljl20 05) find that the abundance profile becomes more 
"boxy" when the diffusion coefficient declines as a function of ra- 
dius. This can be understood, as the large diffusion coefficient near 
the cluster centre causes metals to be rapidly spread to larger radii 
where they accumulate. For the Perseus cluster, there is no clear ev- 
idence that such a model reproduces the observed abundance dis- 
tribution. In Centaurus, however, the high central abundance re- 
gion is more sharply cut off, indicating that such a model might be 
fa voured. Here, we con sider fiducial models of the form proposed 
bv lRebusco et alJ l2005h : 



D = D n 



(7) 



with r\) being a characteristic scale radius at which the diffusion 
coefficient is £>q, and a controlling how centrally concentrated 
the diffusion coefficient is (larger values of a correspond to more 
centrally concentrated diffusion coefficients). However, such mod- 
els may lead to unphysical values of the diffusion coefficient in 
some areas of the cluster and so we impose the additional con- 
straint that D ^ 0. 1 lc s r where c, is the local sound speed. The nu- 
merical factor is taken from the turbulent diffusion simulations of 
iDennis & Chandrar] J2£)05) that are discussed in Section l5"Tl 

Fig. [6] shows abundance profiles with values of a from -1 to 
3 and a scale radius of 25 kpc. It can be seen that the flat inner 
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Figure 6. Variable diffusion coefficient model at an age of 8 Gyr with D = 
4 x lO^cmV 1 ^", r = 25kpc, a = -1 (dots), a = (solid) a = 1 
(dashed) a = 2 (dot-dashed) and background-subtracted model (thick solid 
line). Models with increasing values of a provide an improved fit to the 
data points. 



profile and sharp curvature around 30 — 40 kpc is better matched 
with increasing values of a. This is consistent with a picture where 
the energy driving the turbulent motion is injected at the cluster 
centre, possibly by the inflation and subsequent buoyant rise of the 
radio bubbles. 
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Figure 5. Simulations with constant diffusion coefficients of 2 x 10 27 cm 2 s 2 x 10 28 cm 2 s 4 x 10 28 cm 2 s 1 and 2 x 10 29 cm 2 s ' at times of 1, 2, 4 and 8 
Gyr (dotted, dashed, dot-dashed and solid lines respectively) and background-subtracted model (thick solid line). 



5 TURBULENT HEATING OF CLUSTER GAS 

The dissipation of energy contained in turbulent motions of the 
cluster gas may act as a significant heat source in the ICM. This 
process can provide the kind of spatially extended source of heat 
that is needed to balance cooling over the entire cooling radius 
of the cluster. There are two distinct sources for such energy; the 
merging of sub-clusters and energy injection through the central 
AGN. These processes differ in the expected length scales of the 
energy injection. In the case of mergers, this scale can be compara- 
ble to the size of the cluster, whilst in the case of AGN activity, the 
scale will be much smaller, limited to the size of the radio bubbles. 



Turbulence may also produce heating by mixing of hotter gas into 
the cooler regions, i.e . turbulent conduction (e.g. ICho et alj|2003t 
IVoigt & Fabianl2004) . 



5.1 Heating by dissipation 

On dimensional grounds, a diffusion coefficient may be written as 
D = Ci vl where v and / are velocity and length scales relevant to the 
diffusion process and C\ is a constant, typically of order unity. By 
a similar argument, the simplest possible expression for the rate of 
heating due to dissipation turbulent motions is T — Cipv^ /I, where 
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p is the gas density and v and I are the energy injection scales of 
the turbulence, assumed to be the same as the scales of the diffu- 
sive motions. If we adopt these expressions and assume the cluster 
cooling to be entirely balanced by turbulent dissipation, we may 
write: 

mr)Ar)Hr?= C ^^, (8) 
1(f) 

where A(T,Z) is the cooling functi on. To fix the constants C i 
and C2 we take the values given in iDennis & Chandranl 120051) : 
C\ =0.11 and Ci = 0.4. In practise, Equation|8|is an upper limit 
on the required turbulent heating, as any process that induces tur- 
bulent motion will also provide other heating. Observations show 
both weak shocks and sound waves associated with rising bubbles 
which will provide additional distributed heating. 

Fig. Q shows the calculated length and velocity scales of the 
turbulent motions required to balance cooling, given the variation 
in D with radius. The cooling functi on was calculated using the 
MEK AL plasma emission code fMewe et alJll985l iLiedahl et alJ 
Il995h . assuming that the ratio of all heavy element abundances to 
solar follow that of iron. At the scale radius, the models have a 
scale length of about 4kpc and a velocity of about 150kms _1 . The 
models with a larger than show an initial rise in the length scale 
of turbulent motions — this is over the region where the diffusion 
coefficient is limited by the requirement that D ^ 0.1 lrc, — and 
is the region in which the diffusion approximation is least reliable. 
Outside this inner region, the rapid drop-off in both length and ve- 
locity scales is again consistent with a picture of turbulent energy 
being injected into the cluster core and dissipating away from this 
region. 

Using the calculated length and velocity scales, we may make 
a comparison of the implied turbulent viscosity with the standard 
ISpitzeJ 1 19561) viscosity for a hydrogen plasma, assuming that the 
appropriate turbulent viscosity is of the form = apvl. Using vl = 
D/0. 1 1 and the cluster properties at a radius of 25 kpc, the Spitzer 
viscosity is 2.5 x 10 2 gcm _1 s _1 , and the turbulent viscosity jl — 
4.5 x 10 3 a, so the turbulent viscosity exceeds the ion viscosity for 
any a > 0.06. 

It is notable that, in all the models, the length scale of the 
turbulence rises a factor of more than 3 above the mean ra- 
dius measurem e nt for the Centaurus bubbles of 1.37 kpc from 
iDunn & Fabianl l2004h . If we consider this bubble radius as the 
upper bound on the energy injection radius / max = 1.37kpc, then 
v(r) ^ £>(r)/0.11/ max . Fig.|8|shows the heating as a fraction of the 
cooling rate and the velocity as a fraction of the local sound speed 
assuming / = / max everywhere. The effect on the calculated heating 
rate and velocity distribution is shown in Fig. [8] Since the heating 
rate at fixed D varies as Z~ 4 , fixing the injection scale of the tur- 
bulence to be small causes large increases in the calculated heat 
injected. For models with a > 0, it also requires that the gas ve- 
locity be strongly supersonic in order to maintain the diffusion co- 
efficient in the cluster centre. In practise, the supersonic velocities 
are in a regime where the simple model of turbulence we have used 
does not apply and shock waves would presumably damp such mo- 
tions quickly. Nevertheless, such supersonic motions and large heat 
excesses are not observed in the cluster core, which suggests that 
bubble-injected turbulent motions may not be able to reproduce the 
observed abundance distribution. 

For comparison, Fig. [9] shows th e same information f or the 
Perseus cluster, based on the results of lRebusco et alJ i2005l) . who 
did not find any reason to favour the a > results over the sim- 



ple a = uniform diffusion coefficient. The value of / max is 
fixed at 8.5 2 kpc based on the average inner bubble radius from 
IDunn & Fabianl |2004). In the case of Perseus, the model does not 
produce supersonic velocities, at least for the case of a uniform dif- 
fusion coefficient. Application of the technique to a larger sample 
of clusters would be needed to show whether the problems seen in 
Centaurus are a special feature of that cluster. 

It is clear that reproducing the observed abundance distribu- 
tion with a diffusion model requires both a large diffusion coef- 
ficient at the centre and a relatively small coefficient (D < 4 x 
10 28 cm 2 s~') at radii of about 25kpc to prevent spreading of the 
sharp abundance edge. On the assumption that turbulence is pro- 
duced by the rising bubbles, this situation is hard to reconcile with 
the presently observed bubble properties. This may be because the 
presently observed bubbles are not representative of those produced 
during the lifetime of the cluster; if previous generations of bub- 
bles have been substantially larger, the diffusion coefficient can be 
maintained with lower gas velocities. However the power required 
to inflate the bubbles scales as bubble so tne mcrease m bub- 
ble radius required to bring the diffusive velocity under the sound 
speed under the assumptions used in Fig.|8|would require more than 
an order of magnitude increase in jet luminosity in the models with 
a >0. 

The difficulties may also be avoided if the turbulence is pro- 
duced by a different mechani sm, e.g. focusi ng of sound waves pro- 
duced in sub-cluster mergers iPrinale 1989). Centa urus is believed 
to be undergoing a merger 1 Churazov et alJll999» . If this merger 
is affecting the abundance distribution, for example by causing the 
sharp truncation of a naturally extended profile, comparison of the 
diffusive models to the observed distribution will be misleading. 

If the merger has not affected the abundance distribution, we 
must consider the possibility that some process other than turbu- 
lence has been r esponsible for the observed abundance distribution. 
ICrawford et alJ l2005t) observed spatially coincident Ha and soft 
X-ray filaments surrounding the BCG and conclude that the struc- 
ture of the filaments is inconsistent with turbulence at scales above 
a few kpc. Similar filaments are observed in the core of the Perseus 
cluster iConselice et all200lh and spectroscopic measurements of 
their line-of-sight velocity (Hatch et. al. 2005, submitted) indicates 
that the velocity gradients of the filaments are smooth, which is 
a strong indication that the surrounding medium is not turbulent. 
Velocity measurements of a horseshoe-shaped filament located be- 
low a rising bubble are consistent with laminar flow induced by 
the bub ble rising in a viscous fluid, as proposed by Fabi an et alJ 
(2003b) based on the filament shape. These observations suggest 
filaments may be dragged out of the cluster centre by the rising 
bubbles, a process that has also been implicated in reproducing ob- 
served features in the Perseus abundance map ISanders et all 2005). 
Since these motions are largely radial, they should be more ef- 
ficient at redistributing the metals than turbulent motion and we 
can estimate a diffusion coefficient D ~ Vhnbhipfhnhhie- Usi ng the 
radii and buoyancy velocities from IDunn & Fabianl 12004) as an 
estimate of n^Me an ^ v bubble> we nn d a diffusion coefficient of 
D~2x 10 29 cm 2 r', which will be reduced by the fact that bub- 
ble related flows operate in only a fraction of the cluster core at a 
time. Spectroscopic measurements of the filaments in Centaurus are 
needed to confirm the flow pattern is similar to that in the Perseus 
cluster. 

A second possibility is suggested by Fig.|4] It is clear that the 
the largest discrepancy between the observed abundance distribu- 
tion and the naive abundance-following-light model occurs in the 
highly enriched region inside about 20 kpc. The metal injection into 
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Figure 7. Length and velocity scales needed for turbulent diffusion to balance heating with in a model with Do = 4 x 10 28 cm 2 s rg = 2 5 kpc and a = 
(solid line), a = 1 (dashed line) and a = 2 (dot-dashed line). The cross indicates the position and average radius of the Centaurus bubbles iDunn & Fabian 
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Figure 8. Heating rates as a fraction of the required heating and the velocity as a fraction of local sound speed, assuming D = 0.1 lv/ max (see text for definitions) 
in a model with Do = 4 x 10 28 cm 2 s~', ro = 25 kpc and a = (solid line), a = 1 (dashed line) and a = 2 (dot-dashed line). 



these regions would enhance the cooling rate in these regions, po- 
tentially causing the gas to cool out of the X-ray band. The mass of 
this gas is of the order 10 8 Mq, therefore this would not represent a 
substantial cooling rate over the lifetime of the cluster. 



5.2 Turbulent Conduction 

If we imagine the cluster as being divided into thin shells with 
each of the gas properties being constant within the shell, the 
conductivit y required for heating to balance cooling is given by 
IVoigt & Fabiar] §004 ) : 



These results suggest the diffusion model may be systemati- 
cally overestimating the level of turbulence. It should, however, be 
emphasised that the models of turbulent heating and diffusion are 
very coarse and so these results should be regarded as very approx- 
imate. In particular, the result that turbulent motion on the scale 
size of the bubble requires velocities close to or above the sound 
speed to reproduce the observed iron peak should be regarded as 
demonstrating the limitations of the purely stochastic model of 
metal spreading. 



£/ =1 A(7},Z,-)/i?AV i 



4nr 



(9) 



In the case of an unmagnetised plasma, the theoretical electron con- 
ductivity is given by the Spitzer form: 



""Spitzer = 640 



\m e J 



ke*(kT)i 
e A ZlnA' 



(10) 
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Radius / kpc Radius / kpc 

Figure 9. As Fig.|8]but for the Perseus cluster using the results of lRebusco et alj 12005"): Do = 2 x 10 29 cm 2 s~', ro = 50kpc and a = (solid line), a = 1 
(dashed line) and a = 2 (dot-dashed line). 



where A (not to be confused with the cooling function) is the 
coulomb logarithm: 



lnA«5 37.8 + ln 



IQ S K 



(11) 



In practise the conductivity is likely to be suppressed by the pres- 
ence of magnetic fields which bind electrons to field lines and so 
greatly reduce the effectiveness of transport perpendicular to the 
local field direction. The overall degree of suppression is therefore 
strongly dependent on the field geo metry in the cluster centre. In 
the case of a turbulent magnetic field, Naravan & Medvedev 1 2001) 
find that the effective conductivity is about one fifth the Spitzer 
value. 

In the case of a turbulent g as, the thermal conductivity K is 
given by K" tur b = n e kD Ch o et alj Fig. llOl shows the thermal 

conductivity required to balance conduction in Centaurus with the 
Spitzer conductivity and the turbulent conductivity derived from 
the diffusion model in Fig. [6] The turbulent conductivity is lower 
than the required conductivity everywhere in the cluster. However, 
it appears that Spitzer electron conductivity may be significant in 
the outer regions of the cluster, although this depends on the degree 
of suppression due to the magnetic field. 



6 CONCLUSIONS 

The observed central abundance peak in the Centaurus cluster 
can be approximately reproduced by diffusive spreading of metals 
ejected from the BCG. The boxy abundance profile favours a model 
where the diffusion coefficient declines away from the centre; this 
is qualitatively consistent with the gas motions resulting from the 
rising radio bubbles observed in the cluster core. In particular, re- 
producing the sharp observed abundance drop requires a diffusion 
coefficient of D < 4 x 10 28 cm 2 s~' for r > 25kpc and hence the 
velocity field must be relatively quiet in this region. This is a robust 
result and our main conclusion. 

We have calculated approximate length and velocity scales for 
turbulent motion to heat the cluster. They predict length scales of 
up to ~ 10 kpc and velocities of up to ~ 400 krns -1 in the inner 



an 



x 0.1 



0.01 




T / 10 7 K 



Figure 10. Level of turbulent conduction compared to the required conduc- 
tion for the cluster (solid blue line) and the Spitzer electron conductivity 
(solid black line). Turbulent diffusion is based on the diffusion model with 
D = 4 x 10 28 cm 2 s~ ' , ro = 25 kpc and a = (solid red line), a = 1 (dashed 
line) and a = 2 (dot-dashed line). In all cases, the level of conductivity due 
to turbulent diffusion is below the conductivity required to balance heating. 



~ 25 kpc. That this length scale is larger than the bubble scale is 
a remaining difficulty as a substantial reduction in the length scale 
of the motion implies a reduction of the diffusion coefficient, pro- 
ducing a poorer fit to the data, or a large increase in the turbulent 
velocity. This suggests that the level of turbulence may not be as 
high as predicted. Spreading of metals in the region r < 25 kpc may 
be due to some other process such as the velocity flow behind ris- 
ing bubbles. More work is needed to quantify how effective such 
a mechanism can be. With the next generation of X-ray spectrom- 
eters, direct measurements of turbulent velocities will be possible 
and will place strong constraints on the model. 

The thermal conductivity arising from turbulent mixing has 
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also been calculated and this has been shown to be insignificant 
except, perhaps, in the very central regions of the cluster. 
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